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1  HE  WIGNER  DISTRIBUTION  FUNC1ION 
WITH  MINIMUM  SPREAD 

INTRODUCTION 

A  number  of  advantageous  features  associated  with  smoothing  a  Wigner 
distribution  function  (WDF)  were  discussed  in  a  recent  report  [1].  At  that 
time,  it  was  shown  that  the  WDF  with  minimum  quadratic  spread,  about  the 
line  t  =  bct  in  the  time-frequency  plane,  was  a  two-dimensional  Gaussian 
function,  when  constraints  of  finite  energy  and  mean-square  duration  were 
imposed  [1,  app.  G].  However,  a  more  appropriate  measure  of  spread  about 
the  origin  in  the  t,f  plane  is  adopted  here  and  minimized,  yielding  a  unique 
waveform  and  corresponding  WDF.  Additionally,  a  reward  measure  for 
concentration  is  shown  to  yield  identically  the  same  optimum  WDF. 

An  additional  property  of  smoothing  two-dimensional  WDFs  was  also 
demonstrated;  namely,  if  two  Gaussian  mountains  are  doubly-convolved  with 
each  other,  the  effective  area  of  the  result  is  greater  than  the  sum  of  the 
two  effective  areas,  unless  the  contours  of  both  WOFs  have  the  same  tilt  and 
ratio  of  major-to -minor  axes  [1,  app.  J].  A  quantitative  investigation  of 
the  effect  of  mismatch  in  these  parameters  on  thp  effective  area  is 
conducted  herein. 

It  is  assumed  that  the  reader  is  familiar  with  the  content  and  approach 
of  the  earlier  report;  accordingly,  this  follow-on  effort  will  be  briefer 
and  will  not  review  the  considerable  history  and  background  of  the  WDF. 
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MINIMUM  QUADRATIC  SPREAD 

It  was  shown  In  [1,  (102)  and  (106)]  that  the  short-term  spectral 
estimate  is  equal  to  the  double  convolution  of  the  WOF  of  the  waveform  s(t) 
being  analyzed  with  the  WDF  of  the  weighting  u(t)  employed.  That  is, 

lSu(t>f)|2  -  )jdtl  exP(-i2wfV  s(V  u*(t  -  = 

c  r  ^ 

=  JJ  dt]  df1  Ws(t1,f1)  Wu(t  -  t1 .f  -  f] )  =  Ws(t.f)  ©  Wu(t,f)  .  (1) 

where  ©  denotes  convolution.  Here, 

r 

Wu(t.f)  =  J  dr  exp(  -i2irfr)  u(t  +  f)  u*(t  -  |)  (2) 

is  the  WDF  of  complex  weighting  u(t);  a  similar  definition  holds  for  WDF 
W$ .  (Generalizations  to  non-Wigner  smoothing  functions  for  Wu  are  given 
in  [1 ,  app.  F ] . ) 

Since  the  WDF  of  waveform  s(t)  has  some  good  energy  localization 
properties  (and  some  deleterious  negative  oscillations),  it  is  desired  that 
the  smearing  in  the  t,f  plane,  implied  by  convolution  (1),  be  minimized. 

That  is,  we  would  like  WOF  W^  of  weighting  u(t)  to  be  as  concentrated  as 
possible  about  the  origin  of  the  t,f  plane.  The  ideal  of  an  impulse, 

<S( t ) <5( f ) ,  is  not  a  legal  WDF,  and  must  he  discarded.  Since  the  left-hand 
side  of  (1)  can  never  be  negative,  we  can  be  assured, by  this  smoothing 
procedure  of  two  WDFs,that  we  will  always  get  a  physically-meaningful 
distribution  in  the  t,f  plane;  that  is,  the  smoothed  distribution  will 
always  be  non-negative  for  all  t,f  and  have  a  volume  equal  to  the  energy  of 
waveform  s(t).  For  example,  see  [1.  (Ill)  et  seq.]. 
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PENALTY  MEASURE  AND  SPREAD 


In  order  to  confine  WDF  Wu  near  the  origin,  we  define  a  penalty 
measure  which  is  zero  at  t.f  =  0,0  and  which  increases  quadratical ly  with  t 
and  f.  Namely,  the  penalty  measure  is 

P(t,f)  =  a^t^  +■  4-rr^b^f ^  -i-  4irctf  ,  a,b,c  real  ,  (3) 

and  the  corresponding  spread  of  the  WDF  Wu  is  defined  as 

I  =  dt  df  Wu(t,f)  P(t,f)  .  (4) 

Contours  of  equal  penalty  In  (3)  are  tilted  ellipses  in  the  t.f  plane;  these 
would  be  selected  upon  observation  of  a  calculated  WDF  of  waveform  s(t) 
in  regions  of  interest,  i.e.,  high  activity. 


Therefore,  real  constants  a,b,c  are  presumed  known.  Define  quantity 


Then,  in  order  that  penalty 

P(t,f)  >  0  for  t.f  *  0,0  , 

it  is  necessary  that 


(5) 


(6) 


0  >  0  .  (7) 

2 

The  property  (6)  was  not  satisfied  by  penalty  function  (f  -  Bct)  in 
[1,  app.  G ] ;  that  function  was  zero  all  along  the  line  f  =  3ct,  allowing 
the  WDF  to  become  impulsive  there. 
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We  also  want  WOF  Wu  in  (1)  and  (4)  to  have  unit  volume,  for  two 
reasons.  First  of  all,  this  will  guarantee  that  the  short-term  spectral 
estimate  on  the  left-hand  side  of  (1)  will  have  a  volume  equal  to  the  signal 
energy,  regardless  of  weighting  u ( t )  employed.  Secondly,  without  this 
volume  constraint,  u(t)  and  Wu  would  collapse  to  zero,  giving  a 
meaningless  spread  value  of  I  =  0  in  (4).  Thus  we  require  that 

1  =  jjdt  df  Wu(t,f)  =  fdt  |u(t)|  2  .  (8) 

Subject  to  this  integral  constraint,  we  want  to  minimize  spread  I  in  (4), 
and  find  the  particular  weighting  u(t)  and  corresponding  optimum  WDF  Wu. 
Notice  that  we  are  imposing  no  constraint  of  positivity  on  W^. 


DERIVATION  OF  SPREAD 


Substitute  (3)  into  (4)  to  get  spread 

I  =  J[dt  df  Wu(t,f)  (a2t2  +  4ir2b2f 2  +  4*ctf)  ,  (9) 

where  WDF  W(j  is  given  in  terms  of  u(t)  according  to  (2).  By  using  the 
results  in  [1,  (G-4)  and  ( G-5 ) ] ,  we  can  express  (9)  solely  in  the  time 
domain  as 


I 


a2  Jdt  t2  |u(t)|  2  +•  b2  Jdt  Ju'(t)|  2  +  2c  Jdt  t  Im{u'(t)  u*(tj}  = 

Jdt  [a2t2  Ju(t)|  2  +  b2  |u '  (t)|  2  +  ict  u(t)  u'*(t)  -  ict  u*(t)  u '  ( t )  ] .  (10) 


For  reasons  to  become  apparent  shortly,  define  complex  constant 


B  = 


_  n  +  ic  . 

u2 


(11) 
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then,  by  using  (5),  we  find 

b2|B|2  =  a2.  (12) 

Now  consider  the  quantity 

1  ( t )  r  b2  |u'(t)  +  Bt  u(t)j  2  - 

=  b2  |u' (t)j  2  +  b2  |b| 2t2  [u ( t ) J  2  +  b2Bt  u(t)  u'*(t)  +  h2B*  t  u*(t)  u' (t)  = 

=  b2  |u'(t)|  2  +  a2t2  ju ( t )|  2  +  (VO*  +  ic)t  u(t)  u'*(t)  + 

+  ({O'-  ic)t  u* ( t )  u'(t)  .  (13) 

Comparison  of  (10)  and  (13)  immediately  reveals  that 

Jdt  T ( t )  =  I  +  {O'  Jdt  t  [u(t)  u'*(t)  +  u*( t )  u 1  ( t )  ]  .  (14) 

We  now  integrate  by  parts,  letting 

U  =  t  u(t),  dV  =  dt  u'*(t)  ,  (15) 

to  find  that 


Jdt  t  u(t)  u'*(t)  =  -  j  dt  [u(t)  +  t  u ' (t)  ]  u*(t)  = 

=  -  jdt|u(t)J2  -  jdt  t  u'(t)  u*(t)  .  (16) 

We  presume  that  u(t)  goes  to  zero  at  t  =  ± «o,  consistent  with  energy 
constraint  (8) . 


When  (16)  is  employed  in  (14),  there  follows 

jdt  T(t)  =  I  -  VTT  Jdt  |u(t)|  2 


(17) 
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Thus  the  desired  expression  for  spread  I  is  given  by  (17)  and  (13)  as 

I  =  b2  Jdt  [u‘(t)  +  Bt  u(t)|  2  +  Jdt  Ju(t)|  2  .  (18) 

This  general  result  holds  for  any  weighting  u(t);  it  is  obviously  positive 
in  all  cases ,  since  Q  >  0. 

OPTIMUM  WEIGHTING 

The  last  term  in  (18)  cannot  be  altered;  it  is  equal  to  as  seen  by 
reference  to  constraint  (8).  Furthermore,  the  minimum  value  for  the 
remaining  term  in  (18)  is  zero  and  is  obtained  for  weighting  u(t)  which 
satisfies  the  differential  equation 

u 1 ( t )  +  6  t  u(t)  =  0  for  all  t  .  (19) 

The  only  solution  to  (19)  is 

uQ(t)  =  A  exp(-  |  Bt^)  for  all  t  ,  (20) 

where  complex  constant  A  is  chosen  for  unit  energy,  and  B  is  given  by  (11). 
That  is,  uQ(t)  has  Gaussian  amplitude-modulation  and  linear  frequency- 
modulation.  The  phase  of  A  is  ambiguous. 

The  resultant  minimum  value  of  spread  I  in  (18)  is  obviously 


where  we  employed  (5).  It  is  always  positive,  as  seen  by  reference  to 
requirements  (6)  and  (7). 
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OPTIMUM  WOP 


The  WOF  corresponding  to  optimum  weighting  (20)  is  obtained  by 
substitution  in  (2),  and  use  of  [1,  (H-17)  and  (PI-18)),  as 


WQ(t,f)  =  2  exp 


2  2  2  2  2 
a  t  +  4tr  b  f  j-  4irCtf 


(22) 


The  area  of  the  contour  ellipse  at  the  1/e  relative  level  is  1/2  in  the  t,f 
plane,  as  expected. 


Observe  that  the  numerator  of  the  exp  in  (22)  is  identically  the 
quadratic  penalty  function  P(t,f)  imposed  in  (3).  That  is,  the  contours  of 
optimum  WDP  (22)  are  identical  to  the  contours  of  equal  penalty  of  P(t,f)  in 
(3).  This  result  is  intuitively  satisfying:  the  optimum  WOF  packs  as  much 
volume  inside  a  given  penalty  contour  as  possible,  to  the  extent  that  the 
resultant  WDP  values  are  equal  all  along  that  given  penalty  contour. 


Observe  also,  that  although  positivity  of  the  WDF  W^  was  not  imposed 
as  a  constraint  in  the  minimization  of  spread  I  in  (4)  or  (9),  the  resultant 
optimum  WOP  in  (22)  is,  in  fact,  everywhere  positive.  Although  the  optimum 
weighting  (20)  has  an  ambiguous  phase,  the  optimum  WDP  has  no  ambiguity; 
there  is  a  unique  optimum  WDF,  namely  (22). 
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ALTERNATIVE  REWARD  MEASURE 

Instead  of  penalizing  the  spread  of  WDF  W^  about  the  origin  in  t,f 
space,  we  could  alternatively  utilize  a  measure  which  rewards  concentration 
about  t,f  =  0,0.  In  particular,  consider  reward  function 

R(t,f )  =  exp[-a2t2  -  4ir2b2f2  -  4*ctf]  (22a) 

and  reward  value 

V  =  dt  df  R(t,f )  W  (t,f )  (22b) 

for  WDF  Wu-  The  origin  value  of  R ( t , f )  is  1;  in  order  for  R(t,f)  to  decay 
to  zero  as  t  and/or  f  tend  to  infinity,  we  must  have  condition  (7)  satisfied 
again.  Notice  that  the  contours  of  equal  reward  are  ellipses  in  the  t.f 
plane. 

The  maximization  of  reward  value  V,  subject  to  volume  constraint  (8)  on 
,  is  conducted  in  appendix  A.  It  is  shown  there  that  the  optimim 
weighting  is  again  (20),  and  that  the  optimum  WDF  is  (22).  The  maximum 
value  of  reward  V  is 

Vmax  =  TTfrr  =  rrylt5”!2’-  (22c) 

More  general  results,  for  arbitrary  reward  functions  R(t,f)  in  (22b),  are 
presented  in  appendix  A. 


8 


TR  8317 


GENERALIZATION  TO  SMOOTHED  WDF 


A  general  class  of  distributions*  has  been  presented  in  [2,  (1.7)  and 
(1.8)].  In  current  notation,  that  class  is  given  by  [1,  ( F -1 ) ]  as 

tf 

D(t,f ) =  Wu(t,f)  ®  V2(t,f)  = 

=  du  dT  exp( i2irut  -  i2irft)  ^(v.r)  Q2(v,t)  ,  (23) 

where  WDF  W^  is  given  by  (2),  and  V  (t,f)  is  a  general  two-dimensional 
smoothing  function.  The  complex  ambiguity  function  of  u(t)  is 

Xy(».T)  =  Jdt  exp(  -i2iri>t)  u(t  +  f)  u* ( t  -  |)  ,  (24) 

while 

q2(u,r)  =  JTdt  df  exPM2™t  +  12*fl)  V  (t.f)  (25) 

is  a  double  Fourier  transform  of  the  smoothing  function  V?.  Observe  that 
if  there  is  no  smoothing,  then 

V2(t,f)  =  6(t)  6(f) 
q2(»,r)  =  1  for  all  u,T 

D(t.f)  =  Wu(t,f)  .  (26) 


*This  section  is  based  upon  a  suggestion  by  Leon  Cohen,  Hunter  College,  New 
York,  NY,  that  the  optimum  WDF  results  here  actually  apply  to  a  wider  class 
of  distributions. 
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Now  it  is  shown  in  appendix  B  that  the  following  second  moments  of 
generalized  smoothing  distribution  D  can  be  expressed  in  terms  of 
derivatives  of  ")C  and  q,,  at  the  origin: 

ffdt  df  t2  D(t,f)  =  -  ~  [%JU(0,0)  q2(0,0)  +  2^(0, 0)  qj(0.0)  + 

J  4* 

+  ^(0.0)  qpO.O)]  , 

(fdt  df  t  f  0(t,f)  =  ~  ftJ’fO.O)  q2(0,0)  e}^(0,0)  q”(0,0)  + 

4ir 

+  xj(0.0)  qj(0,0)  +XU(0.0)  q2T(0t0)]  , 

J[dt  df  f2  0(t,f)  =  -  ~  (%JT(0,0)  q2(0,0)  +  2X^(0, 0)  qj(0,0)  + 

+  XU(0,0)  qjr(0,0)]  .  (27) 

Here,  for  example,  superscript  v  denotes  a  partial  derivative  with  respect 
to  v,  which  is  then  evaluated  at  the  origin  v,T  =  0,0. 

If  follows  immediately  that  if  origin  value 

q2(0.0)  =  1  ,  (28) 

and  if  the  five  origin  derivatives 

q"(0,0)  =  qj(0,0)  =  qjw(0,0)  =  qjr(0,0)  -  q2T(0,0)  =  0  ,  (29) 

then  (27)  reduces  to  the  moments  that  would  have  resulted  from  employing  the 
no-smoothing  result  (26)  in  (27).  Thus,  distributions  D(t,f)  resulting  from 
(23),  with  properties  (28)  and  (29)  for  q^,  have  the  same  second  moments 
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as  the  WDF  W  (t,f).  Hence,  the  spread  Ip  of  distribution  D(t,f)  is 
given  by  (see  (9)) 

IQ  =  [Jdt  df  Oft.fi  (a2t2  +  4*Vf2  *  4irctf)  = 

=  JJdt  df  Wu(t,f)  (a2t2  +  4ir2b2f2  +  4»ctf)  =  I  ,  (30) 

which  is  exactly  the  spread  I  of  WDF  Wu(t,f).  That  is,  smoothed 
distribution  D(t,f)  in  (23)  has  the  same  spread  as  WDF  Wu(t,f),  when 
smoothing  function  V2(t,f)  (actually  transform  q^)  satisfies  the 
properties  in  (28)  and  (29).  Notice  that  these  properties  are  considerably 
less  restrictive  than  requiring 

P2(».0)  =  Q2(0,T)  =  1  for  all  v.X  ,  (31) 

which  arises  when  one  is  interested  in  maintaining  the  marginals  [2,  (1.6)]. 

We  must  also  observe  from  (23)  that  the  volume  under  generalized 
smoothing  distribution  D  is  equal  to  the  product  of  the  volume  under  W^ 
and  the  volume  under  V^.  But  the  latter  quantity  is  unity,  by  virtue  of 
(28).  What  all  this  means  is,  that  if  we  minimize  spread  Ip  in  (30), 
subject  to  a  unit  volume  constraint  on  D,  the  end  result  is  precisely  (18) 
and  (20),  and  the  optimum  WDF  W^  is  again  given  by  (22).  The 
corresponding  distribution  D  is  obtained  by  substitution  of  (22)  into  (23) 
and  specification  of  the  complete  or  q^  functions.  The  properties  in 
(28)  and  (29)  are  not  sufficient  to  completely  specify  q^  or  D;  all  that 
is  specified  by  (28)  and  (29)  are  the  second  order  moments  of  D  in  (27).  It 
should  also  be  noted  that  all  the  conditions  in  (29)  cannot  be  met  by  the 
general  tilted  Gaussian  q?  function  employed  in  [1,  (F-9)  and  sequel  to 
(F— 12)]. 
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SENSITIVITY  TO  MISMATCH 

In  [1,  app.  J],  it  was  shown  that  if  two  Gaussian  mountains  are 
doubly-convolved  in  x,y  space,  the  effective  area  of  the  resultant  is 
greater  than  the  sum  of  the  individual  areas,  except  when  the  two  elliptical 
contours  have  the  same  tilt  and  the  same  ratio  of  major-to-minor  axis  (shape 
factor).  Here,  we  wish  to  investigate,  quantitatively,  the  increase  in 
effective  area  above  the  minimum  value,  when  the  tilt  and  shape  factors  are 
not  at  their  optimum  values.  This  situation  can  arise  when  observation  of 

WDF  of  waveform  s(t)  is  contaminated,  in  a  particular  region  of 
interest  in  the  t,f  plane,  by  interference  effects,  thereby  making 
estimation  of  the  tilt  and  the  shape  factor  of  the  elliptical  contours 
somewhat  inaccurate. 

The  general  situation  is  considered  mathematically  in  appendix  C. 

Ellipse  1  has 

area  A1 ,  tilt  0^  shape  factor  F1  ,  (32) 


while  ellipse  2  has 

area  A?,  tilt  0?,  shape  factor  F?  . 


(33) 


The  ratio 


A,  +  A. 


(34) 


is  presented  in  (C-13)  in  terms  of  a  number  of  auxiliary  quantities. 
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The  initial  example  we  consider  is  where  ellipse  1  has  seven  different 
areas,  namely 

A1  =  .5, 1,2, 3, 4.5, 6  01  -  J  F1  =  2  •  (35) 

Ihe  tilt  is  fixed  at  w/4  radians  and  the  shape  factor  at  2.  On  the  other 
hand,  ellipse  2  has 


A2  =  2  02  =  -  J  to  J  F2  =  2  .  (36) 

That  is.  the  shape  factor  is  perfect  at  F?  =  F  =  2,  but  the  tilt  is 
swept  over  a  u/2  range  (greater  discrepancies  than  w/2  lead  to  obvious 
periodicities  and  symmetries  centered  about  S2  =  0^  as  well  as  about 
02  =  0^  ±  */2  and  about  0?  =  0^  +  *) .  The  situation  under 
investigation  is  depicted  in  figure  1,  where  ellipse  2  is  dotted. 

The  effect  of  mismatch  in  tilt  is  presented  quantitatively  in  figure  2. 

As  expected,  ratio  (34)  is  1  at  0=0^*  ir/4,  regardless  of  area  A] . 

Ihe  most  degradation  (upper-most  curve)  is  realized  for  A1  =  2,  i.e.,  when 
the  areas  of  the  two  ellipses  are  equal.  The  maximum  increase  in  area  is 
only  25  percent,  when  02  is  off  by  ir/2  radians;  however,  if  the  shape 
factor  is  significantly  larger  than  1,  the  sensitivity  to  the  tilt  would  be 
much  greater,  as  figure  1  shows. 

The  final  example  utilizes  the  exact  same  parameter  values  as  (35)  for 
ellipse  1,  while  ellipse  2  has 


F2  =  2  to  6. 


(37) 
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Now  the  tilt  is  perfect  at  0^=8^=  »/ 4,  but  the  shape  factor  F? 
varies  above  the  best  value  of  2.  The  situation  Is  depicted  in  figure  3, 
where  ellipse  2  is  again  dotted. 

Ratio  (34)  Is  plotted  in  figure  4  versus  the  shape  factor  F  .  Again, 
the  upper-most  curve  corresponds  to  the  case  where  A^  -  A^  =  2.  There 
is  no  need  to  compute  ratio  (34)  for  <  F^  =  2,  because  the  values  for 
F2  =  F^r  are  the  same  as  those  for  F2  =  F^/r.  Additional  cases  of 
interest  can  be  investigated  by  use  of  the  program  listed  in  appendix  C. 
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Figure  2.  Area  Ratio  (34)  for  Fj  =  2,  F2  =  A?  =  2,  Pi  =  w/4 
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SUMMARY 

The  most  compact  WDF  Wu  that  can  be  used  for  two-dimensional  smoothing 

of  a  measured  WDF  W  is  a  Gaussian  function  in  two  variables,  when  the 
s 

measure  of  spread  is  quadratic  in  the  time  and  frequency  variables  t  and  f, 
or  the  reward  measure  is  exponential.  Furthermore,  this  two-dimensional 
convolution  guarantees  a  non-negative  modified  distribution,  since  the 
result  is  equivalent  to  a  short-term  spectral  estimate.  Extensions  to  a 
particular  class  of  generalized  distriDutions  yields  the  same  optimum  WOF. 
The  corresponding  waveform  has  Gaussian  amplitude  modulation  and  linear 
frequency-modulation. 

The  additional  smearing  caused  by  mismatched  smoothing  functions  to  the 
true  parameters  of  a  measured  WDF  has  been  investigated  numerically  for  a 
few  examples,  and  found  not  to  be  overly  sensitive  to  the  exact  values. 
However,  the  multitude  of  parameters  has  prevented  simplification  of  the 
area  spread  factor;  accordingly,  a  program  allowing  calculation  of 
particular  cases  is  included  to  allow  for  further  investigation. 

The  WDFs  for  the  Hermite  functions  of  order  n  are  given  in  closed  form, 
in  terms  of  a  Laguerre  polynomial  of  order  n.  This  result  is  extended  to 
cross-WDFs  in  appendix  A;  in  this  manner,  we  can  investigate  the  WDF  of  an 
arbitrary  waveform  when  expanded  in  a  weighted  sum  of  Hermite  functions, 
including  linear  frequency-modulation. 
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APPENOIX  A.  MAXIMIZATION  OF  REWARD  VALUE 


We  want  to  find  that  WDF,  W  (t,f),  which  is  maximally  concentrated 
about  the  origin  in  t,f  space,  where  the  measure  of  reward  for  concentration 
is 

R(t,f )  =  exp[-a2t2  -  4ir2b2f 2  -  4*ctf]  ,  a.b.c  real  .  (A-l) 

Thus,  the  maximum  reward  occurs  at  the  origin, 

R(0.0)  =  1  ,  ( A -2 ) 

and  the  contours  of  equal  reward  are  ellipses  in  the  t,f  plane.  In  order 
for  R ( t , f )  to  decay  to  zero  as  t  and/or  f  tend  to  infinity,  we  must  have 


0  >  0  , 


( A -3 ) 


where 


n  2,2  2 

0  =  a  b  -  c 


( A  -4 ) 


The  reward  value  associated  with  WDF  W  is  the  real  quantity 

u 


V  =  JJ  dt  df  R(t.f)  Wu(t.f)  , 


( A  -  5 ) 


which  we  wish  to  maximize,  where 


Vt.f)  =  [  dr  exp( -i2irfr)  u{t  +  u*(t  -  |) 


( A -6 ) 


in  terms  of  weighting  u(t).  We  must  constrain  the  volume  of  Wu>  in  order 
that  V  in  (A-5)  not  tend  to  infinity  as  u(t)  is  simply  increased  in  level. 
Thus,  we  have  integral  constraint 


1  =  JJdt  df  Wu(t,f)  =  Jdt  | u ( t ) | 


( A- 7 ) 
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ALTERNATIVE  FORM  FOR  V 

If  we  substitute  (A-6)  in  (A-5),  there  follows 

V  =  ft  dt  dt  r(t,T)  u(t  +  f)  u*(t  -  f)  ,  ( A-8) 

where 

r(t.t)  =  Jdf  exp( -i2irf t)  R(t,f)  .  (A-9) 

A  more  useful  alternative  form  for  (A-8)  is 

V  =  j^dx  dy  K( x , y )  u(x)  u*(y)  ,  (A-10) 

where  kernel 

K(x,y)  =  r ,  x-y^  .  (A-ll) 

EIGENFUNCTIONS  OF  K 

In  this  and  the  following  subsection,  kernel  K  is  Hermitian,  but 
otherwise  arbitrary;  it  is  not  limited  to  form  (A-ll)  with  (A-9)  and  (A-l). 
Suppose  and  £0^  are  the  eigenvalues  and  eigenfunctions  of  kernel 
K;  i .e. , 

jdx  K(x.y)  en(x)  =  Xn  0n(y)  for  n  =  0,1,2,  ...  ,  (A-12) 

where  xQ  >  >  x^  ...  ,  and 

f dx  0*(x)  0  (x)  =  6  (A-13) 

J  n  nr  nm 
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Then  the  kernel  can  be  expanded  according  to 

oO 

K(x,y)  =  Xn  0*(x)  «n(y)  .  (A-14) 

n=0 

Also,  there  follows  immediately 

JJdx  dy  K(x,y)  an(x)  0*(y)  =  .  (A-15) 

EXPANSION  OF  u 

Suppose  we  expand  weighting  u  in  a  series  of  eigenfunctions  of  Hermitian 
kernel  K: 

u(x)  =  ^  gn  0n(x)  ,  where  gn  =  Jdx  u(x)  a*(x)  .  (A-16) 

n=0 

Then  general  reward  expression  V  in  (A-10)  becomes 

co 

V  =  |  dy  u*(y)  Jdx  K(x,y)  9n  an(*)  = 

n=0 

aC  OQ 

■  5  9„  [dy  «*(y)  k„  »„(y)  -  ^  |9„|2  ,  (A-.7) 

n=0  n=0 

where  we  used  (A-l?)  and  (A-16).  At  the  same  time,  the  energy  of  u  in 
(A-16)  is 

E,  =  [dx  }u(x)|  2  =  ]>  |gn|2  .  (A-18) 

n=0 
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Now  if  the  energy  E^  of  u  is  constrained  at  1 ,  as  in  (A-7),  then  the 

best  choice  of  coefficients  K\  to  maximize  V  in  ('*-17)  is,  since 

X  >  x  >  X„  . . . ,  obviously 
o~l-2 

|gQ|  =  1  and  gn  =  0  for  n  >  1  .  (A-19) 

That  is,  the  optimum  weighting  is 

u0(xi  =  0o(x)  exp(ie)  ,  (A-20) 

where  constant  e  is  arbit-ary,  while  the  maximum  reward  is 


V  =  x  . 
max  o 


( A  -  2 1 ) 


lhat  is,  the  zero-th  order  eigenvalue  and  eigenfunction  of  general  Hermitian 

kernel  K  in  (A  12)  are  the  solutions  to  the  problem  of  interest  here,  namely 

maximization  of  reward  value  V  in  (A-10)  by  choice  of  weighting  u.  For  a 

general  kernel,  a  recursive  numerical  procedure  could  be  employed  on  (A-12) 

to  determine  \  and  0  ,  if  desired, 
o  o 


The  formulation  in  these  last  two  subsections  is  actually  general  enough 
to  cover  the  earlier  penalty  function  considered  in  (3)  et  seq.  The  only 
difference  is  that  the  eigenvalues  {x^  now  i nc rease  with  n,  and  we  must 
select  the  eigenfunction  corresponding  to  the  minimum  eigenvalue,  in  order 
to  realize  the  least  penalty.  This  approach  is  the  subject  of  appendix  D. 
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SPECIAL  CASE  OF  EXPONENTIAL  REWARD 


We  now  specialize  the  general  results  of  the  previous  two  subsections  to 
the  reward  function  (A-l).  Substitution  in  (A-9)  and  use  of  (A-4)  yields 


r(t’x)  =  exp 


?  1  ? 

Qt  +  ~-T  -  ictT 


( A -22) 


compare  [1,  (F-9)  and  (F-12)].  Then  (A-11)  immediately  gives  Hermitian 
kernel 


K(x.y)  -  afc  expf- 


2y?,b 


4b 


( A -23 ) 


where 


D  =  Q  +  1  +  i  2c  . 


(A-24) 


At  this  point,  we  refer  to  Mehler's  expansion  [3,  (67)]  to  obtain 
(after  some  labor) 

oO 

K(x,y)  =  xn  «*(x)  dp(y)  , 
n=0 


( A — 2  5 ) 


where 


and 


,  (i  -VoY 
n  (i  +  Vo1) 


n+1  • 


«n(x)  =  A  exp 


IaI2  =  ^ 
‘A>  yrb  • 


(A  -26) 


( A  -27 ) 


(A  -28) 
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The  function  Hen(x)  is  the  Hermite  polynomial  [4,  22.2.15].  It  is  easily 
verified  that  (A-27)  satisfies  orthonormality  relation  (A-13). 


OPTIMUM  WEIGHTING 


Since  Q  >  0  by  ( A -3 ) ,  the  eigenvalues  in  (A-26)  satisfy 

x  >  >  x„  . . .  .  Therefore,  the  maximum  reward  is 

o  1  2 

■  ttw  ■  Try ' 

and  the  corresponding  weighting  is 


u0(t)  *  «0U)  »  A  exp 


l 


J?  fQ*  +  ic 
2  b2  . 


( A— 29 ) 


( A— 30 ) 


from  (A  27)  and  (A-28).  This  is  identical  to  (20)  combined  with  (11). 
Therefore  the  optimum  WDF  is  again  (22)  for  reward  measure  (A-l),  as  well  as 
penalty  measure  (3).  The  waveform  in  ( A -30 )  has  Gaussian  amplitude 
modulation  and  linear  frequency-modulation. 


HIGHER -ORDER  HERMITE  FUNCTIONS 


For  n  >  0,  the  reward  values  fxn\  in  (A-26)  are  all  less  than  optimum 

value  xq.  We  have  succeeded  in  obtaining  these  explicit  values  without 

having  to  evaluate  the  WOFs  of  the  corresponding  Hermite  waveforms  in 

(A-27).  We  now  rectify  this  situation.  The  WDF  of  t  (t)  in  (A-27)  is 

n 

given  by  integral 
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Wn(t,f)  -  Jdx  exp(-12*fT)  #n(t  +  |)  0*(t  -  |)  = 

J^r  Iac  *«!>[-«*«  ■  H - 1)2]* 
*  H%((1  *  ?))  H*„(F£  -  ?))• 


( A -31 ) 


where 


B  ,  M±JC  F  =  2l^  |A|2-^ 

b2  b  ’  ‘A1  "  yrb  • 


(A-3 2) 


Now  a  more  general  integral  result  already  exists  in  closed  form;  from 
[5,  p.  292,  (30)],  we  have,  in  a  form  more  useful  for  present  purposes, 

Jdx  exp( -  2  x?  +  dX)  Hem(b  +  x)  Hen(b  -  x)  = 


=  (-l)m  m!  (b  -  a)n  m  l^n  m^(b^  -  a^)  exp(a*V2)  for  m  <  n  , 


( A  -33) 


where  Lma^x^  is  the  generalized  Laguerre  polynomial  [4,  22.2.12].  When 

( A— 33 )  is  used  on  ( A— 32 ) ,  there  follows,  for  the  WOF  of  waveform  0  (t)  in 

n' 

(A-27),  the  compact  result 


Wn(t,f)  =  (-1)  2  Ln(2U)  exp(-U)  , 


( A -34 ] 


where 


2+2  *  2,2,2  . 

at  *■  4*  b  f  +  4irctf 


( A-35) 


This  result  reduces  to  (22)  for  n  =  0.  Again,  contours  of  equal  values  of 
the  WOF  are  ellipses  in  the  t,f  plane. 
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CROSS-WDFs 

Suppose  a  general  waveform  u(t)  is  expanded  in  a  set  of  orthonormal 
Hermite  functions  with  linear  frequency-modulation  (a  positive  real,  B  real) 

1/4  .  , 

dn(t)  =  (J)  exp[-  j  (a  +  iB)tN  Hen(yT7  t)/{nT ,  (A-36) 

according  to 

to 

u(t)  =  un  *n(t)  .  ( A— 3 7 ) 

n=0 

Then  the  WDF  of  u(t)  becomes 

Wu(t,f)  =  Jdr  exp( -i2»fr)  u(t  +  J)  u*(t  -  f)  = 

<x3 

=  u  u*  W  (t,f)  ,  (A-38) 

cL—  m  n  mnv  ’  '  ’  v  ' 

m,n=0 

where  cross-WDF 

Wmn(t,n  =  [dr  exp(-i2*fT)  dm(t  +  f)  i*(t  -  J)  .  ( A— 39 ) 

When  (A-36)  is  substituted  in  (A-39),  and  (A-33)  is  utilized,  the 
cross-WDF  can  be  expressed  as 

Wmn(t.f>  ■  L<n-'">(|z|2)  «p(-|z|2/2)  for  .  <  n  . 

( A-40) 

where 

z  =^[a t  +  i(2wf  +  Bt)]  , 

\z\2  =  l  [(a2  +  02)t2  +  4»2f2  +  4irBtf]  .  ( A -41 ) 
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These  results  generalize  [6,  pp.  456-7]  and  [7,  p.  547]. 

The  origin  value  of  (A-40)  is 

%(°-°)  =  2<-1>m  4mn  •  <A~42> 

consistent  with  unit  energy  of  ( t )  and  their  even  or  odd  character.  The 
cross-WDF  in  (A-40)  is  a  function  only  of  the  three  variables  m,n,z,  where  z 
is  the  complex  combination  in  (A-41).  The  parameters  a  and  B  in  (A-36)  are 
perfectly  general;  when  they  are  specialized  to  match  (A-27),  and  when  we 
set  m  =  n,  then  (A-40)  reduces  to  (A-34).  Equations  (A-38)  and  (A-40) 
afford  a  direct  calculation  of  the  WOF  of  a  general  waveform  u(t),  once  the 
coefficients  are  determined  by 

Un  =  I dt  u(t)  0n(t)  •  (A-43) 
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APPENDIX  B.  MOMENTS  OF  DISTRIBUTION  D 

The  generalized  smoothing  distribution  D  is  given  by  (23)  in  terms  of  a 
double  Fourier  transform  of  product 

P(v,T)  =  X  (v.D  q2(v.r)  .  (B-l) 

Therefore,  the  inverse  relation  is 

P( v.-ir)  =  JJ  dt  df  exp( -i2irvt  +  12*fT)  D(t,f)  .  (B-2) 

If  we  let  superscript  v  denote  a  partial  derivative  with  respect  to  v,  there 
immediately  follows  from  (B-2), 

P(0,0)  =  JJ  dt  df  D(t,f ) 

Pu(0,0)  =  -i2ir  JJ  dt  df  t  D(t,f ) 

PT(0,0)  =  i2ir  JJ  dt  df  f  D(t,f ) 
pvv(0,0)  =  -u2  JJ  dt  df  tz  D(t,f ) 

PVT(0,0)  =  4ir2  JJ  dt  df  t  f  D(t,f ) 

P^fO.O)  =  -4ir2  JJ  dt  df  f2  D(t,f )  .  ( B-3) 

When  these  relations  are  written  out  explicitly  in  terms  of  and 
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q2,  according  to  (B-l),  we  find  that  the  moments  of  D  are 

JJ  dt  df  O(t.f)  =  >£(0,0)  q2(0,0) 

JJ  dt  df  t  0(t,f)  =  ^  Dcj<0.0)  q2(0.0)  +*u(0.0)  qj(0,0)] 

\\  dt  df  f  0(t,f)  =  ^  [xj(0,0)  q2(0.0)  +  >^(0.0)  q’J(O.O)] 

JJ* dt  df  t?  D(t,f )  =  -  -*7  EjcJ*(°.0)  q2(0.0)  +  2X^(0. 0)  q”(0,0)  + 

4 w 

+  *u(0,0)  q^CO.O) ] 

f f  dt  df  t  f  O(t.f)  =  ~  (*JX(0.0)  q2(0,0)  *  *J<0,0)  qj(0,0)  t 
J  4ir 

+  X”(0.0)  qJ(O.O)  +%u(0,0)  qJT(0,0)] 

jjdt  df  f2  D(t,f )  =  -  ~~  CXurr(0,0)  q2(0,0)  +  2*^(0, 0)  q£(0.0)  + 

4ir 

+  ZU(0,0)  qXT(0,0)  ]  .  (B-4) 

Since  q^  is  a  double  Fourier  transform  of  V^,  of  exactly  the  same  form 
as  (B-2),  it  follows  immediately,  by  similarity  to  (B-3),  that  the  required 
derivatives  of  q^  in  (B-4)  can  be  found  from  smoothing  function  as 
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q2(°.°)  =  Jjdt  df  V2(t.f) 

q”(0,0)  =  — i 2ir  J{  dt  df  t  V2(t,f) 

q{(0,0)  =  i2ir  JJ  dt  df  f  V2(t,f) 

q2v(0,0)  =  -4*2JJ  dt  df  t2  V2(t,f) 

q^O.O)  =  4*2{Jdt  df  t  f  V2(t,f) 

q^O.O)  =  -4ir2Jjdt  df  f2  V2(t,f)  .  (B-5) 

COMPLEX  AMBIGUITY  FUNCTION  PROPERTIES 

The  required  derivatives  of  ^  in  (B-4)  can  be  determined  from 
definition  (24).  We  list  them  here  for  completeness  and  future  reference: 

Xu(0,0)  =  jdt  |u(t)|2 
Xj(O.O)  =  -i 2ir  {dt  t  Ju( t) j  2 
Xj(0.O)  =  i  {  dt  Im  [u 1  (t)  u*(t)J 
*”U(0,0)  =  -4*2  Jdt  t2ju(t)J  2 
^(O.O)  =  2«  jdt  t  Imfu'(t)  u*(tj] 

X^O.O)  =  -  Jdt  Ju '  (t)|  2  .  (B-6) 
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These  quantities  are  all  real,  with  the  exception  of  the  two  single 
derivatives,  both  of  which  are  purely  imaginary.  These  second-order 
derivative  values  of  ^  can  be  expressed  solely  in  terms  of  u(t)  and  u'(t). 

Since  we  can  express  complex  ambiguity  function  ^  in  terms  of  the  WDP 

W  according  to 
u  3 

^(v.T)  =  JJdt  df  exp(  -i2wvt  +  i2irfx)  Wu(t  ,f )  ,  (B  -7) 

it  readily  follows  from  ( 8-6)  that 


Jj  dt  df  Wu(t,f)  =  J  dt)u(t)]  2 
55  dt  df  t  Wu(t,f)  =  Jdt  1 1  u  ( t )  |  2 
JJdt  df  f  W^t.f)  =  ^  Jdt  Imfu'(t)  u*(t^ 

5J  dt  df  t2  Wu(t,f)  =  Jdt  t2  ]u(t)|  2 

5]  dt  df  t  f  Wu(t,f)  -  ^  J  dt  t  Im(u'(t)  u*(t)] 

55  dt  df  f?  Wu(t*f)  =  ^2  5dt  lu’  (t)J  2  •  ( B  —8 ) 
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APPENDIX  C.  GENERAL  TILTED  ELLIPSE 

It  will  be  convenient  to  be  able  to  specify  the  area,  tilt,  and  shape 
factor  of  an  ellipse  directly,  instead  of  trying  to  solve  for  these 
quantities  from  the  general  form 

12  12  r 

j  ax  +  ~  by  +  \faF  pxy  =  1  (C-l) 

employed  in  [1,  (J-2)].  Accordingly,  as  done  In  [1,  app.  0],  we  employ  the 
rotated  coordinates  depicted  in  figure  C-l  below.  The  equation  of  the 
el  1  ipse  in  x ' ,y '  space  is 


Figure  C-l.  Rotated  Coordinate  Axes 
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But  since  the  area  of  this  ellipse  is 


while  its  shape  factor  is 


A  =  ir  x'  y'  , 
o  o 


X' 

F  =  — 

y '  ’ 
*o 


( C  -  3 ) 


(C-4) 


it  is  a  simple  matter  to  find  that 


A  F  =  ir  x^2  ,  A/F  =•  tr  y^2  , 


( C  —  5 ) 


leading  to  the  desirable  form 


(C-6) 


Furthermore,  the  coordinate  axes  in  figure  C-l  are  related  according  to 


x'  =  xC 


y1  =  -xS 


+  ys"| 

+  yc  J 


C  =  cos { B )  ,  S  =  s i n ( B )  . 


(C-7) 


Substitution  in  (C-6)  yields 


-  x2^-  +  FS2)  +  l  y2(j^  +  FC2)  xy  SC  Q  -  f)  = 


(C-8) 


which  is  of  the  form  (C-l)  under  identifications 


|i,FS3),  „.&(£*  FC2)  , 

fTT5  wil#  ’  ■  sc(f  - F)  • 


with  y  =  SC | 


(C  9) 


Once  area  A,  tilt  8,  and  shape  factor  F  are  specified,  (C-9)  affords  a  ready 
calculation  of  a,b,p;  quantities  C  and  S  are  given  by  (C-7).  Since 
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p  =■•  $in(e)  In  [1,  (J-6)],  we  have 


sin(e)  = 


* 


*■  y 


=5p,  cos(e) 


(C-10) 


which  are  needed  below. 

In  order  to  distinguish  the  two  Gaussian  mountains  being  doubly 
convolved  in  [1,  (J-2)],  we  label  them  with  subscripts  1  and  2, 
respectively,  thereby  obtaining 


a "  *r(h  *  f,s>)  • b  *  *  Firii  • s’  ■  s,n,v  •  c>  ■  tos,“') 


Y,  -  S 


!C1  |F~  -  F,\  .  sin(e)  J..x~yr,  cos(e)  -  T 

'  1 VI  7  Vi  +  y‘  Vi  + 


(C-ll) 


and 


/c2  \  (s2 

C  =  +  F2S2j  ,  d  =  “1^  +  F2C2J  ,  S2  =  sin(B2)  ,  C?  =  cos(02) 


y2  -  s2c2 


( h  '  Fj)  ■ 


sin(0)  -  n—  ~2r<  cos(0)  = 


fT  VI  +  y 


r  • 
2 


(C-l  2) 


We  are  now  in  a  position  to  evaluate  the  effective  area  of  the 
resultant  convolution;  namely  from  [1,  ( J-9) -( J-l 1 ) ] ,  we  have 


JH 


A1  K  A2  ^b^cosfe)  +  Vcd’cos(^) 


(C-l 3) 
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where 

2  2  r 

0  =  ab  cos  (e)  +  cd  cos  (0)  +  ad  +  be  -  2 yabc d'  sin(e)  sin(a)  .  (C-14) 

The  minimum  value  of  (C-13)  is  1,  attained  when  shape  factors  f^  = 

and  tilts  B^  =  B0.  More  generally,  when  we  specify 

,B^ , for  ellipse  1  , 

A2*°2,F2  f0r  e11ipse  2  ’  (C-15) 

equations  ( C -1 1 )  and  ( C -1 2 )  allow  for  evaluation  of  all  the  parameters 
needed  In  (C-13)  and  (C-14).  A  sample  program  in  BASIC  is  attached. 
Subroutine  E  computes  a,  b,  sin(e),  cos(ft)  as  given  by  (C-9)  and  (C-10)  in 
terms  of  given  area  A,  shape  factor  F,  and  tilt  B  (=B). 
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APPENDIX  D.  KERNEL  APPROACH  TO  PENALTY  FUNCTION 

The  general  formulation  in  (A-11)  through  (A-21)  will  be  applied  in 
this  appendix  to  the  penalty  function  (3): 

P(t,f)  =  a2t2  +  4ir2b2f 2  4-  4*ctf  .  (D  1) 

Substitution  in  (A-9)  (in  place  of  reward  R)  yields 
r(t,T)  =  Jdf  exp( -i2trfr)  P(t,f)  = 

=  Jdf  exp(  -i2irfT)  (a2t2  +  4ff2b2f2  +  4»ctf)  = 

=  a2t24(t)  -  b24"(X)  +  i2ct4' (T)  .  ( D-2) 

Then  kernel  K  follows  from  (A-11)  as 

K(x,y)  =  ,  x  -  y)  = 

2 

=  (X  +  y)2  4(  x  -  y)  -  b24"(x  -  y)  +  ic  (x  +  y)  4’(x  -  y)  ,  (D-3) 

which  is  Hermitian. 

The  integral  equation  (A-12),  that  must  be  solved,  can  be  simplified  by 
use  of  the  facts  that 

4  a2  |  dx(x  +  y)2  4(x  -  y)  «n(x)  =  a2y2ein(y)  , 

-b2  [  dx  4"(x  -  y)  en(x)  =  -b2«T(y)  , 
ic  J  dx  (x  +  y)  4 1 ( x  -  y)  »n( x)  =  -ic[2ye^(y)  +  0n(y)]  .  (D-4) 
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where  the  last  two  results  are  obtained  by  integration  by  parts.  Then 
( A-l 2 )  yields  differential  equation 


b2«r(y)  +  i2cy  0^{y)  +  (\n  +  ic  -  a2y2)  *n(y)  =  0 


(D-5) 


If  we  try  solution 


«0(V)  =  A  exp  ^  3/ 


(0-6) 


in  (0-5),  we  find  it  to  be  acceptable  if  we  take 


•  \-ff- 


(0-7) 


These  results  agree  with  (11)  and  (20),  as  expected.  To  find  the  general 
solution  of  (D-5),  we  try  solution  form 


O(y)  =  exp(~  ^  By2)  H( y)  , 


(0-8) 


with  B  still  given  by  (0-7).  This  form  in  (D-8)  is  no  loss  of  generality 
since  H  is  still  arbitrary.  Use  of  (0-8)  in  (D-5)  results  in 


b 2 H " ( y )  +  2y  H'(y)  ( -b2B  +  ic)  + 


+  H(y)  ( -b2B  +  b2B2y2  -  i2cBy2  +  x  +  ic  -  a2y2)  =  0 


(0-9) 


When  the  value  for  B  in  (0-7)  is  utilized,  (0-9)  simplifies  to 


b  H»(y)  -  2TPJ*  y  H'(y)  +■  (x  -  )  H(y)  =  0  . 

(As  a  partial  check,  if  H(y)  =  A,  then  \  =  Vcf,  as  in  (0-7).) 


(0-10) 


Now,  in  (0-9),  let 


H(y)  --  G(  Fy)  ,  H '  (y)  =  F  G'(Fy)  ,  H"(y)  -  r  G "  ( F  y )  ,  (0-11) 
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where  F  is  arbitrary  for  the  moment,  thereby  obtaining 


b2F2  G"(Fy) 


Now  let  x  =  Fy  to  get 


6"  ( x) 


-  2V0’y  F  G'(Fy)  +  (X  -  tfT)  G(Fy)  =  0 


-  ^  x  G '  ( x )  +  —-p’ctx)  =  0  . 


b  F 


b  F 


If  we  now  let  (without  loss  of  generality) 

r  21/2Q1/4 


then  (0-13)  simplifies  further  to 

G"(x)  -  x  G 1  ( x )  +  ^^GIx)  =  0 
We  now  appeal  to  [4,  22.6.21]  and  observe  that  if 

X-^JGT. 


2  fir 


=  n  =  integer  , 


then  a  solution  of  (D-15)  is 


Also,  (0-11)  yields 


while  (0-8)  gives 


with 


G(x)  =  Hen(x)  ,  xn  =  VTJ»  (1  +  2n) 


H(y)  =  G(Fy)  =  Hen(Fy)  , 


«n(y)  =  A  exp(-  j  By2)  Hen(Fy)/Vn?  , 


V7T  +  ir  -l/2nl/4  -  nl/4 

r  =  IS  1  jc  F  =  2  _  Q_ i.(2  _  Q  , 

h2  ’  b  •  ’ 


(0-12) 


(0-13) 


(0-14) 


(0-15) 


(0-16) 


(0-17) 


( 0—18) 


(0-19) 


(0-20) 
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where  the  unit  energy  norma  1 i zation  of  «n  has  been  imposed.  The 
corresponding  eigenvalue  follows  from  ( D -1 7 )  as 

Xn  =  VTn  +  2")  =\/a2b2  -  c2' (1  +  2n)  .  (0-21) 

The  minimum  obviously  occurs  for  n  =  0. 

Result  (0-19)  agrees  with  (A-27).  However,  the  xn  given  here  by 
(0-21)  differs  from  that  given  by  (A-26),  because  we  are  solving  for  the 
minimum  penalty  here  versus  the  maximum  reward  there. 
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